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EFFECT  OF  EVAPORATION-CONDENSATION 
ON  PHOTOACOUSTICS  OF  AEROSOLS 

FINAL  REPORT 

SUMMARY 

This  preliminary  investigation  of  relative  humidity  effects  on  light  absorption  by  aerosols  has 
been  very  successful.  The  Desert  Research  Institute  (DRI)  took  part  in  the  Southern  Great  Plains 
Radiation  Measurement  and  took  simultaneous  measurements  of  photoacoustic  light  absorption  and 
light  absorption  measured  by  filter  methods  as  a  function  of  relative  humidity.  The  variation  in 
relative  humidity  was  controlled  artificially.  This  work  is  described  in  detail  in  Section  I. 

A  detailed  theory  of  the  effect  of  evaporation-condensation  on  photoacoustic  output  has  been 
developed  at  the  University  of  Mississippi.  This  theory  is  valid  over  a  wide  span  of  particle  sizes, 
particle  densities  and  operating  frequencies.  This  work  is  presented  in  Section  EL 

Richard  Raspet,  Pat  Amott  and  Hans  Moosmtiller  met  at  DRI  in  May  2000  to  coordinate  the 
research  effects  and  to  plan  future  cooperative  research.  The  measurements  and  theory  will  be 
combined  in  the  near  future  in  the  preparation  of  papers  for  publication,  of  experimental  design  and  of 
proposals  for  future  research. 

1.  RELATIVE  HUMIDITY  DEPENDENCE  OF  LIGHT  ABSORPTION 
The  photoacoustic  instrument  is  used  to  obtain  light  absorption  by  aerosols,  with  applications 
relevant  to  radiation  transfer  in  the  atmosphere,  or  more  specifically,  electro-optical  imaging  and 
sensing.  The  photoacoustic  instrument  has  the  potential  for  obtaining  the  relative  humidity 
dependence  of  light  absorption,  whereas  more  common  filter  methods  cannot  provide  a  quantitative 
measure  due  to  moisture  deposition  on  the  substrate.  For  dry  aerosols  (RH<50%),  the  sound  generated 
in  the  instrument  is  due  to  the  conversion  to  heat  of  the  light  absorbed  by  the  aerosols.  The  aerosols 

are  sufficiently  small  that  this  heat  transfers  entirely  to  the  surrounding  air  during  the  acoustic  period. 
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However,  for  moist  aerosols,  vapor  may  be  released  from  the  heated  aerosols,  in  addition  to  heat.  In 
the  general  case,  the  sound  pressure  level  is  due  to  both  heat  and  mass  transfer. 

The  contribution  of  mass  transfer  to  the  sound  pressure  level  measured  in  a  photoacoustic 
instrument  was  investigated  by  performing  measurements  and  theoretical  analysis.  The  photoacoustic 
instrument  was  transported  to  the  Department  of  Energy,  Southern  Great  Plains  Atmospheric  Radiation 
Measurement  in  north  central  Oklahoma  (abbreviated  as  SGP).  An  “aerosol  trailer,”  among  others,  is 
located  at  the  SGP,  and  the  purpose  of  this  trailer  is  to  measure  optical  and  physical  properties  of 
atmospheric  aerosol.  In  particular,  personnel  from  NOAA  in  Boulder,  CO  operate  a  system  to  vary  the 
relative  humidity  (RH)  of  the  incoming  air  from  40%  to  95%.  This  RH  conditioner  has  been  evaluated 
by  NOAA,  and  provided  a  very  convenient  means  to  accomplish  the  experimental  component  of  the 
effort.  Ancillary  instruments  also  present  in  the  trailer  were  useful  in  quantifying  aerosol  properties 
before  the  humidity  conditioner.  The  relative  humidity  dependence  of  light  absorption  has  not  been 
directly  measured  in  the  past  due  to  the  use  of  filter-based  methods.  The  photoacoustic  instrument  was 
described  previously,1,2  and  the  calibration  can  be  verified  with  nitrogen  dioxide  gas.3 

The  data  set  from  this  project  is  outstanding.  DRI  personnel  were  present  for  the  project  (the 
entire  month  of  March,  2000),  and  the  aerosol  trailer  was  operated  during  this  time  to  ensure  good  data 
quality.  A  synopsis  follows. 


1.  Aerosol  light  absorption  measurements  with  3  instruments. 

a.  Normal  psap  filter  based  instrument  used  at  SGP. 

b.  Extra  psap  sitting  on  the  humidified  inlet. 

c.  DRI  photoacoustic  instrument  on  the  humidified  inlet. 

2.  Light  absorption  for  diy  and  humidified  aerosols  was  measured.  Light  absorption 
measurements  are  among  the  more  difficult  to  do.  A  great  combination  of  aerosols  was 
available: 

a.  Cases  of  very  clean  air  (background) 

b.  Cases  of  long-range  transport  of  probably  vehicle  exhaust  (tractors,  interstate,  roads,  etc, 
mostly  diesel  perhaps?) 
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c.  Cases  of  smoke  plumes  from  farmer's  fires. 


3.  Findings  to  date  are: 

a.  The  photoacoustic  data  correlates  well  with  the  DRY  psap  measurements.  However,  the 
exact  value  for  light  absorption  is  not  settled  yet.  (Slope  of  correlation  is  not  one,  but 
appears  to  be  around  1 .6,  with  the  photoacoustic  instrument  reporting  less  light  absorption 
than  the  other  instrument). 

b.  The  photoacoustic  data  show  no  increase  of  light  absorption  with  increase  of  RH,  and 
perhaps  a  decrease.  (Theory  suggests  perhaps  a  factor  of  2  increase  for  carbon  cores 
coated  with  an  aqueous  layer.  But  are  light  absorbing  aerosols  hygroscopic?) 

c.  The  humidified  PSAP  shows  erroneously  high  light  absorption  as  the  RH  increases  due  to 
the  filter  taking  on  aqueous  solution  from  scattering  aeroso.  When  the  RH  goes  back 
down,  this  psap  shoots  erroneously  negative  as  well.  It  is  very  good  that  the  SGP  crew 
uses  a  dry  psap  for  their  measurements  on  a  routine  basis.  It  also  calls  into  question  all 
measurements  ever  taken  with  a  psap  at  elevated  RH,  including  those  done  on  aircraft. 

To  evaluate  the  light  absorption  dependence,  the  ratio  of  light  absorption  obtained  by  the 
photoacoustic  instrument  to  that  of  the  dry  psap  was  computed  and  averaged  over  the  entire  data  set. 
This  ratio  is  deemed  the  f(RH)  function  for  light  absorption,  and  is  a  dimensionless  value.  If  f(RH) 
increases  as  the  RH  increases,  than  the  conclusion  is  that  light  absorption  increases  with  RH  as  the 
aerosols  swell  and  moisten  with  increased  RH.  It  may  also  be  that  the  chain  aggregate  structure  of  the 
carbonaceous  aerosol  collapses  upon  humidification.  Another  potential  effect  is  that  RH  increases 
cause  aerosol  losses  in  the  sample  lines  to  increase.  It  should  be  kept  in  mind  that  the  analogous  f(RH) 
for  scattering  by  aerosols  is  in  the  range  of  2  to  4  for  the  RH  around  95%. 

Figure  1  shows  the  measured  f(RH)  for  a  1  fxm  inlet  on  the  sample  line.  First,  the  average  ratio  is 
around  0.6,  indicating  that  the  psap  light  absorption  is  larger  than  the  photoacoustic  value.  It  would  be 
unity  if  both  instruments  reported  the  same  light  absorption.  This  calibration  issue  is  likely  to  be 
traced  to  the  psap  calibration  method.  Next,  note  that  the  f(RH)  function  decreases  with  RH.  This 
decrease  could  be  due  to  at  least  the  following  three  factors: 


1.  The  sample  lines  begin  to  absorb  particles  as  the  RH  increases. 
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2.  Mass  transfer  is  less  efficient  at  producing  sound  than  is  heat  transfer,  so  the  photoacoustic 
instrument  could  report  less  light  absorption  for  elevated  RH.  Theory  should  provide  a 
bound  for  this  efficiency  difference. 

3.  Aerosol  morphological  changes  could  lower  the  mass  absorption  efficiency  if  the  chain 
aggregate  particles  of  soot  collapse  upon  humidification. 


Figure  2  shows  the  measured  f(RH)  for  a  10  pm  inlet  used  on  the  sample  line,  allowing  larger 
particles  through  the  system.  Note  that  the  f(RH)  function  drops  off  faster  with  increasing  RH  for  this 
inlet  (i.e.  larger  aerosol)  than  it  does  for  the  1  pm  inlet. 


RH  (%) 

Figure  1.  Ratio  of  photoacoustic  and  dry  psap  light 
absorption  as  a  function  of  relative  humidity.  The  inlet 
line  had  a  1  pm  cut  to  prevent  larger  particles  from 
entering.  The  line  is  a  least-squares  fit  to  the  data,  and  is 
included  to  illustrate  the  trend  (not  to  hypothesize  a 
linear  relationship.) 


RH  (%) 

Figure  2.  Same  as  Fig.  1,  though  the  inlet  size  cut 
allowed  particles  smaller  than  10  pm  to  pass  through  the 
instruments.  Note  the  relatively  smaller  ratio  in  this 
figure  at  high  RH,  when  compared  with  that  of  Fig.  1. 


In  summary,  the  relative  humidity  dependence  of  light  absorption  was  measured  as  a  function  of 
relative  humidity  using  a  photoacoustic  instrument,  and  for  dry  aerosol  using  a  filter-based  instrument. 
The  data  suggests  that  RH  increases  are  associated  with  a  decrease  in  light  absorption.  Three 
hypotheses  have  been  provided  for  the  mechanism  of  decrease.  The  data  set  will  be  further  analyzed, 
and  compared  with  theory. 
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2.  THEORY  OF  THE  EFFECT  OF  EVPORATION-CONDENSATION 
ON  PHOTOACOUSTIC  OUTPUT 


2.1  Introduction 

DRI  has  developed  photoacoustic  spectrometers  for  the  measurement  of  light  absorption  by 
atmospheric  aerosols.1’2  A  laser  source  is  chopped  near  the  resonant  frequency  of  an  acoustic  resonator 
containing  the  sample.  The  light  absorbed  by  the  aerosol  in  the  beam  heats  the  surrounding  gas,  driving 
the  acoustic  wave  in  the  resonator.  The  sound  pressure  level  in  the  resonator  is  a  measure  of  the  light 
absorption.  The  results  of  this  technique  have  been  compared  to  those  from  an  aethalometer,  a  device 
that  measures  the  absorption  of  particles  gathered  on  a  filter  medium.  DRI  is  interested  in  the  effects 
of  relative  humidity  on  the  adsorption  of  light  by  aerosols  and  on  the  effect  of  volatile  components  on 
the  calibration  of  the  photoacoustic  instrument. 

M.B.  Baker  analyzed  the  effect  of  evaporation-condensation  on  the  photoacoustics  of  volatile 
aerosols  in  a  non-resonant  cell.4  Baker’s  analysis  assumes  that  the  aerosol  particles  do  not  interact  and 
that  the  photoacoustic  cell  is  non-resonant.  In  this  paper,  we  allow  for  high  aerosol  particle  number 
density  effects  and  perform  the  Green’s  function  source  analysis  as  appropriate  for  resonant  cells.  The 
analysis  is  based  on  the  low  reduced  frequency  method  used  previously  to  evaluate  the  effects  of 
evaporation-condensation  on  sound  propagation  in  wet  walled  tubes.5,6 

2.2  Theory 

The  DRI  photoacoustic  cell  is  a  flow  through  cell.  The  laser  beam  is  much  smaller  than  the  area 

of  the  resonate  tube.  We  therefore  assume  that  the  residence  time  of  the  droplets  in  the  beam  is  small 

compared  with  the  time  for  significant  DC  mass  evaporation  to  occur.  The  laser  excitation  can  be 

considered  as  a  DC  component  producing  steady  evaporation  and  heating  and  an  AC  component  that 

will  alternately  heat  and  cool  the  droplet.  Our  analysis  will  concentrate  on  the  AC  component.  In 

addition,  we  assume  that  each  droplet  of  radius  r0  is  surrounded  by  similar  neighbors  and  that  drop 
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inertia  is  high  enough  that  each  droplet  supplies  heat  and  vapor  into  its  neighborhood  of  radius 
a  =  l/«1,/3  where  n  is  the  aerosol  particle  density  in  m‘3. 

The  governing  equations  are  the  Navier-Stokes  equation,  the  continuity  equation  for  each 
component  and  for  the  mixture,  the  equation  of  state  for  the  mixture,  the  entropy  equation  for  the 
mixture,  and  the  diffusion  equation  for  the  mixture.  In  linearized  form  these  are:7'8 


P„^-=-Vp+AtV2v+fj 0+ipWv), 
St  3  J 


^+P?(V-v,)=  0, 

^+P;(v.v2)=o, 

J^  +  p„(Vv)  =  0, 

~gl+Y~[nkYt~’^2T  +  kTnkTy(vl-v2)=0, 


—  =  nk  —  +kT 
dt  dt  ° 


1  ^Pl.+.  1  dPl 


{mx  dt  m2  dt  j 


and 


vi-v2  =D,2 


-mx)vP~— — -^vr 

P2  Pi  PoPo  nxn2  Ta 


(1) 

(2a) 

(2b) 

(2c) 

(3) 

(4) 


(5) 


In  these  equations,  variables,  pltp2’  the  density  of  the  gas  and  of  the  vapor;  p=  p2  +  p2  is  the 
density  of  the  mixture;  Vj,v2  are  the  hydrodynamic  velocity  of  the  gas  and  of  the  vapor; 
v  =  (pjV,  +  p2v2  )/p  is  the  hydrodynamic  velocity  of  the  mixture;  p  is  the  pressure  of  the  mixture;  and 
T  is  the  temperature  of  the  mixture.  In  the  above  equations,  constant  |i,  is  the  viscosity  of  the  mixture; 
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(3  is  the  bulk  viscosity  of  the  mixture;  m1/m2  are  the  molecular  mass  of  gas  and  vapor;  n1/n2  are  the 
number  density  of  the  gas  and  of  the  vapor;  n  =  n1  +  n2  is  the  total  number  density;  are  the  ratio 
of  specific  heats  of  the  gas  and  of  the  vapor;  ny/(y-l)=  nlyl/(yl  -l)+n272/(72  is  the  ratio  of 
ncp/R  for  the  mixture;  cp  is  the  molar  heat  capacity  at  constant  pressure;  R  is  the  universal  gas 

constant;  K  is  the  thermal  conductivity  of  the  mixture;  D12  is  the  mass  diffusion  coefficient;  kT  is  the 
thermal  diffusion  ratio,  as  defined  in  Ref.  8,  page  541;  and  k  is  the  Boltzmann’s  constant. 

The  linearized  heat  equation  for  the  droplet  is  given  by, 

^ipLCpLTL)-KLV2TL=Q(l),  (6) 

where  pL  is  the  liquid  density;  cpl  is  the  liquid  heat  capacity  at  constant  pressure;  Kl  is  the  heat 
conductivity;  Tl  is  the  temperature  of  the  liquid;  and  Q(t)  is  the  energy  absorbed  from  the  laser  per  unit 
time  per  unit  volume  in  the  drop.  Assuming  single  frequency  excitation,  we  set Q(t)  =  Q0e~,ax ,  where 
Qa  =  a  1 0  with  the  absorption  coefficient  for  the  drop  given  by  a  and  the  intensity  of  the  laser  light 

given  by  f,.  We  assume  that  the  absorption  of  light  energy  is  small  enough  so  that  the  drop  is  evenly 
heated.  We  assume  also  that  the  resulting  diffusion  and  heat  flow  will  be  predominately  radial  and 
define  two  dimensionless  variables, 


r 

a 


the  dimensionless  radius,  and  the  reduced  frequency. 


c 


(7) 


(8) 
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Assuming  that  all  the  variables  have  an  e'1C0t  time  dependence,  where  to  is  the  angular  frequency  of  the 
modulated  laser  beam,  we  introduce  normalized  variables  as  follows: 


p=P°c2  Lo+p*e-i0*), 

(9a) 
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T  =  Tol^.0+T*e~i<a), 

(9b) 

TL=T0(l  +  T*e-iw), 

(9c) 

v  =  c  v*e',ox , 

(9d) 

and 

v,—  v2  =  cVV®,  (9e) 


where  c  is  the  sound  speed  in  the  mixture. 

Using  the  definitions  above,  the  Navier-Stokes  equation  becomes 


y  dii  T]2di i[t]2  Bt] 


(10) 


where  AM  =  aJ°vj  P +— M  •  Note  that  this  definition  of  the  dimensionless  viscous  wave  number  is  not 


the  same  as  the  shear  wave  number  defined  in  Ref  5.  The  entropy  equation  is  given  by. 


r  4  n2  dm  dji  r  n 


pcoc 

where  hr  =  a, - — .  The  equation  of  state  is  combined  with  the  continuity  equations  of  each 

V  k 


component  to  form  the  equation  of  state  in  terms  of  p\  T*,  v*  and  V* : 
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The  final  equation  results  from  combining  the  diffusion  equation  with  the  continuity  relations; 


_*  /  d  f  1  d(/72v*V)  Q  n  /  \dp*  nzkT  dT* 

— - - - — >  — - — — [m,  - w,  I— - - - =0, 

^DdrlW  8j7  J  tfblPo  8rl  n\ni  _ 


where  XD  =  a  Jyp  •  The  temperature  equation  in  the  droplet  is. 


r  -±±±\  n2  ?£\  =  iQ° 

1  Al  T]1  dTj{n  dri  )  cpLpLG)T0 


(17) 


ndT 

■  Q - h 

dri 


i£lV  + 


n2y 


J  +  0 

V  4, 


_3_ 

3  TJ 


i  dfov) 

(tj2  3tj 


:«l”2 

«Po 


(m,  -m2)— 
077 


_l4?v) 

tj2  dr] 


=  0. 


At  this  point,  we  note  that  the  reduced  frequency,  Q,  is  small,  even  for  low  number  density  and 
high  operating  frequency.  (Q  =  0.04  for  n  =  350  particles/cm3  and  f  =  2000  Hz).  Accordingly,  we 
drop  all  terms  in  Q  from  the  equations: 


n 

r*_,-_L_L± 

f  2 dr*) 

«l«2 

y-i 

/l£  T]2  077 

dr] 

\  / 

r/2  dr] 

,np0 

(mi  -m2)-kT 


(18) 


Q.  n2kr  dT*  ( 

v *  1  a 

i  din2v')) 

X2D  dr] 

X2d  drj 

y  J_ 

(19) 


and 


n3r*  ,  a 

077  07? 


i  afev)) 

772  077 


-mj- 


np„ 


df] 


1  d(n2v') 

r l2  ft) 


=  0. 


(20) 


The  term  containing  v*  can  be  eliminated  by  combining  Eqs.  (18)  and  (20): 


dT *  i  d 

_1 _ 3_ 

f „ 2  dr*  j 

'  ,  .y-l*r  d 

i  4v)’i 

drj  a2  an 

rj2  dr] 

r  ndr] 

W  *  J 

(21) 


Equations  (19)  and  (21)  are  coupled,  second  order  differential  equations  for  - —  and  V*.  The 

3tj 

coupling  is  provided  by  the  thermal  diffusion  ratio,  kr,  that  appears  in  the  entropy  equation  from 
diffusion  across  the  control  volume  and  in  the  diffusion  equation  for  temperature  driven  diffusion. 

These  equations,  combined  with  boundary  conditions  and  the  temperature  equation  in  the  drop, 
represent  a  complete  solution.  The  effect  of  the  coupling  on  air-water  mixtures  has  been  investigated 
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in  Ref  6.  For  photoacoustics  in  the  atmosphere,  the  mole  fraction  of  vapor  will  be  very  low  (less  than 
0.05)  and  the  effect  of  the  coupling  of  the  diffusion  and  thermal  wave  numbers  low.  In  Ref  6,  the 
calculated  fractional  change  in  thermal  and  diffusive  wave  numbers  for  water  vapor  in  air  is  less  than 
1CT4  for  temperatures  up  to  30°C.  The  use  of  the  low  reduced  frequency  approximation  and  the  modal 
analysis  of  Ref  6  demonstrates  that  Baker’s  use  of  simplified  hydrodynamic  equations  is  justified  for 
this  analysis  for  higher  frequency  devices.  The  boundary  conditions  for  the  solution  of  this  problem 
are  developed  within  these  approximations.  Diffusion  and  thermal  effects  are  assumed  to  dominate  the 
acoustic  terms  in  the  droplet  neighborhood. 


2.3  Boundary  Conditions 

The  boundary  conditions  at  the  neighborhood  boundary  (r  =  a)  are 


V*(1)  =  0, 


and, 


dT 


* 


dr] 


71=1 


=  0. 


(22a) 


(22b) 


At  the  drop  radius,  r0,  the  temperature  must  be  continuous;  the  heat  generated  by  condensation  at 
the  surface  must  be  carried  away  by  the  heat  flux  into  the  gas  and  droplet;  and  the  diffusion  velocity 
can  be  calculated  from  the  temperature  at  the  drop  radius.  We  require  that  the  temperature  in  the  drop, 

T*l  ,  and  the  gas  mixture,  T* ,  be  equal  at  the  drop  radius  (r|  =  r<>/a): 


T 


* 


(23) 


Next,  we  require  heat  flux  balance  at  the  drop  radius.  In  general,  this  can  be  written  as, 
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(24) 


q  =  K\T-KLVTL, 

where  q  is  the  heat  flux  generated  at  the  boundary.  For  a  system  undergoing  a  phase  change  at  the 

wall,  the  heat  flux  generated  is  written,  q  =  tin nux  =  £p2  v2(~)>  where  p°  is  the  ambient  vapor 

a 

r 

density;  v2(— )  is  the  radial  velocity  of  the  vapor  evaluated  at  the  wall;  and  £  is  the  latent  heat  of 
a 

vaporization.  Thus,  our  boundary  condition  for  the  heat  flux  at  T]  =  rQ  / a  becomes, 


£P°2v2=kVT-klVTl.  (25) 

Use  of  the  definition  of  the  radial  difference  velocity,  V* ,  yields  the  following  equation  for  the  heat 
flux. 


(26) 


Temperature  fluctuations  at  the  drop  radius  induce  vapor  pressure  fluctuations.  The  Clausius- 
Clapeyron  equation,9 


P  Vapor  ra/  ~  Po  exP 

*  /  n 


£m. 


T  T, 


REF 


l=Va 


(27) 


relates  the  vapor  pressure  to  the  temperature.  An  expression  for  the  density  of  species  2  near  the  tube 
wall  can  be  derived  from  the  ideal  gas  law  and  Eq.  (27): 


Tref 

T 


exp 


im2 

rl 

1  11 

R 

T 

\ 

Tref 

(28) 


Linearizing  with  T  =  To(l  +  T*e~'6>t)  yields. 
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(29) 


1  ,-.0 
rn  /  —  Py 

f 

1  + 

(im2  ^ 

>1 

T*e-imt 

lr>=/a  H2 

,R7o  , 

> 

where  T0  is  the  ambient  temperature  of  the  system.  The  equation  of  continuity  can  be  used  to  write 
Eq.  (29)  in  terms  of  V*  and  T* : 


j__a_ 

V  dn 


=n\ 

'  Lrri!  j' 

T* 

II 

ft- 

^R70  J 

v=°/ 


(30) 


The  divergence  of  v*  is  dropped  from  this  boundary  condition  since  v*  at  the  boundary  is 
smaller  than  F*by  the  ratio  p2/  P  and  v*  is  principally  an  acoustic  wave  with  small  spatial  variation 
compared  to  diffusion  and  heat  conduction.10  The  fourth  boundary  condition  is  that  the  time  rate  of 
heat  flow  out  of  the  drop  at  r0  must  equal  the  time  rate  of  change  in  internal  energy  in  the  drop 
neighborhood, 


Anr, 


2  k  dT* 


a  8rj 


=  ~i0}i0laCvPT*(TlY™Z112dTl . 


(31) 


Again,  this  boundary  condition  is  only  correct  if  the  acoustic  energy  emission  is  negligible  with  respect 
to  the  internal  energy. 


2.4  Solution 

The  set  of  differential  equations  to  be  solved  are. 


dT *  i  d 

drj  drj 

,72  dtl 

l7  3ti)\ 

(32) 
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tl  =DjoyilLT])+i 


cplPl°>To 


A,  B  and  D  are  constants  to  be  determined  from  the  boundary  conditions.  fx  {jiXrj)  and  /0  (V/Arj) 
are  defined  as, 
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where  r\0  =r0/a.  The  functions  fi  and  fo  are  chosen  so  that  they  approach  one  for  the  small  drop, 
large  neighborhood  limit.  This  facilitates  comparison  with  Baker’s  theory.  Recursion  relations  for  the 
spherical  Bessel  functions  give, 


(V/An))=  i tfriofo  (ViAr?), 

T]  dri 


(40) 


di  n0 


(41) 


and 


d_ 

dll 


=  -iX2fl(j~iXri). 


(42) 


Integration  of  Eq.  (35)  gives, 


T*{ri)  =  ~rioA  fo  (V*V?)+  C , 


(43) 


where  C  is  an  additional  undetermined  constant.  This  term  arises  due  to  the  finite  neighborhood 
assumption.  Using  this  form  in  Eq.  (31)  gives, 


2 

AAKr°K~fx {fiXrTi0 )=  icot  cvp[-rjoAf0 (^7})+  cJfcraV^ . 
a  0 


(44) 


Performing  the  integration  yields, 


A  /1  {fiXprio )  =  -icocv  p 
a 


-^-fxy~iXrrio)+^j-{l-Jll) 


(45) 
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Solving  for  C  we  find, 


C  = 


3  iArc 
cocvpo 4 


(46) 


Usually  Tjl  is  negligible.  C  tends  to  zero  as  a  goes  to  infinity.  The  equation  for  the  temperature  T*  is 
given  by, 

T*(n)  =  -V oA  [/o  ('TM'Tti)-  8  /i(V//V7o  )]>  (47> 

x_Vo_ 

3//t*  T  Y  3  ic  - 

where  8  - - Y - This  term  is  small  for  all  aerosols  of  interest.  - «  4x10  for 

acvp  a3  \-rjl  cvp 

saturated  air/water  at  30°C.  The  spatially  constant  temperature  fluctuations  in  the  neighborhood  of  the 
drop  are  negligible.  In  the  following,  we  set  8  =  0 .  The  temperature  matching  boundary  condition, 
Eq.  (23),  becomes 

DjJ\iaLrlo )+  ,Qo  rr,  =  -M/o (J'Lt’Io )•  (48) 

cplPlc°T0 

The  heat  flux  boundary  condition,  Eq.  (26),  becomes 

“  ^~Bfi  (VUD70)=  —  Afx  (yfiATT}0  )+— Dj\{fiALTj0  ).  (49) 

T0  a  a 

The  boundary  condition  that  expresses  the  vapor  pressure  in  terms  of  the  surface  temperature  of  the 
drop  gives. 
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Eliminating  D  from  Eqs.  (48)  and  (49),  we  find  that 


The  product  A ltj0  is  usually  small;  small  argument  expansions  can  be  used  to  express  Eq.  (51)  in 
simpler  form. 


Let  c 


'  =  s  i  +  —  {al rj0  )2 .  The  second  expression  results  from  the  small 

4iALTj0  JoWM-lVo)  15 


argument  expansion  and  demonstrates  that  c'  can  be  taken  as  1  unless  Altj0  is  of  0(1).  For  100  pm 
drop  of  water  at  2000  Hz  this  value  is  2.3,  at  20  Hz,  Baker’s  assumed  value  it  is  0.07.  Gathering  terms, 


we  find, 


,  .cpLPL6>r02c' 

1  - 

3k 


Qoca 
3  T0 


(52) 


From  Eq.  (50), 


(53) 


Eq.  (52)  is  solved  for  A: 


(54) 


(55) 


J.  * 

These  constants  may  be  combined  with  Eqs.  (35)  and  (36)  to  evaluate  -  and  V  .  For  the  present 

a??  ar* 

investigation  of  photoacoustics,  we  evaluate  the  acoustic  source  terms  rather  than  the  variables - 


and  V 


The  wave  equation  with  mass  injection  and  heat  injection  is  given  by11 


cz  dr 


where 


+^t-h 

dt[  c  )  l  c2  ) 


M  is  the  mass  injected  per  unit  time  and  volume;  H  is  the  heat  energy  added  per  unit  time  and  volume. 
The  mass  injected  per  unit  time  per  unit  volume  is  calculated  by  the  product  of  p2v2  at  the  drop 
surface  times  the  surface  area  divided  by  the  neighborhood  volume  with  v2  =  -cV* .  This  gives 


fo  , 

/i  (V7az??70  ) 

fo 

a  kT„  E  RT 


The  heat  energy  injected  per  unit  time  per  unit  volume  is  determined  by  summing  the  heat  flux  over 
the  surface  divided  by  the  neighborhood  volume: 


„j!k)es 


a  i  a 


For  weakly  absorbing  drops,  Q0  is  the  product  of  the  beam  intensity  and  the  liquid  absorption  in 
inverse  meters.  The  optical  absorption  length  due  to  the  aerosol  is  given  by  the  absorption  in  the 
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droplet  times  the  volume  ratio  ( r0faf  .  The  phase  of  f{r,t)  will  be  the  phase  of  the  resulting 
pressure  wave  with  respect  to  the  phase  of  the  laser  intensity. 

2.5  Analysis 

Equations  (58)  and  (59)  along  with  Eq.  (55)  have  been  written  to  facilitate  comparison  with 
Baker’s  low  frequency,  small  particle,  large  spacing  results.4  Her  results  are  recovered  by  setting 

c'  =  1, 

fo{fi^Trlo  )  fl  (jihp'Ho  )  1 
f\  {Ji^rVo  )  /o  ) 

and 

•  c plPl aro  c  fo{^^rrlo)  q 

3k-  /j  {ji^jTj0  ) 

Eq.  (60a)  expresses  the  condition  that  the  droplet  radius  is  smaller  than  a  thermal  penetration  depth  in 
the  droplet;  that  is,  the  temperature  is  spatially  uniform  in  the  drop.  Eq.  (60b)  expresses  two 
conditions.  The  first  is  that  the  finite  drop  neighborhood  does  not  affect  the  heat  and  mass  transfer  at 
the  drop  surface.  This  condition  will  be  met  if  the  neighborhood  radius  is  more  than  a  thermal 
penetration  depth  in  the  gas.  The  second  condition  inherent  in  Eq.  (60b)  is  that  XLr\0  be  small  enough 
so  that  the  leading  terms  in  the  small  argument  expansion  for  h0  and  hx  are  sufficient.  This  is  due  to 
the  normalization  of  /0  and  fx  retaining  only  the  leading  terms  (see  Eqs.  (38)  and  (39)). 

Figure  3  displays  the  magnitude  and  phase  of  the  ratio  in  Eq.  (60b)  versus  the  neighborhood 
radius  for  two  particle  sizes,  1  |im  and  10  (im  for  typical  frequencies  for  Baker’s  study  and  for  DRI 
instrumentation.  The  ratio  is  significantly  different  from  one  only  when  the  neighborhood  radius  is  the 


(60a) 

(60b) 

(60c) 


2 

size  of  a  thermal  penetration  depth.  At  2000  Hz,  this  only  occurs  for  particle  densities  of  240  cm' 
corresponding  to  a  particle  spacing  of  10'  m. 


Figure  3a.  Magnitude  of  the  ratio  in  Eq.  (60b)  as  a 
function  of  neighborhood  radius  for  different  operating 
frequencies  and  drop  sizes.  2000  Hz  and  10  micron  sized 
drop  (solid  line);  200  Hz  and  0.1  micron  sized  drop 
(dashed  line);  20  Hz  and  1  micron  sized  drop  (dotted  line). 
The  ambient  air  pressure  is  1  atm  and  the  ambient 
temperature  is  290  K. 


Figure  3b.  Phase  of  the  ratio  in  Eq.  (60b)  as  a  function  of 
neighborhood  radius  for  different  operating  frequencies 
and  drop  sizes.  2000  Hz  and  10  micron  sized  drop  (solid 
line);  200  Hz  and  0.1  micron  sized  drop  (dashed  line);  20 
Hz  and  1  micron  sized  drop  (dotted  line).  The  ambient  air 
pressure  is  1  atm  and  the  ambient  temperature  is  290  K. 


The  last  condition,  Eq.  (60c),  is  violated  in  higher  frequency  instruments.  The  magnitude  and 
phase  of  f0{jiATVo)/ fS^^rVo)  f°r  the  same  cases  as  Fig.  3  is  displayed  in  Fig.  4.  These  ratios  are 
near  constant  and  real  for  particle  densities  less  than  240  cm'3.  The  10  pm  drop  at  low  frequencies  is 
affected  by  the  next  term  in  the  expansion  of  h0  and  hx .  The  magnitude  of  this  term  is  0(1)  none-the- 
less.  The  frill  term  in  Eq.  (60c)  can  have  a  significant  effect  on  the  acoustic  level  for  the  high 
frequency  device  and  is  frequency  dependent.  This  term  is  due  to  the  change  in  magnitude  and  phase 
of  the  surface  temperature  due  to  heat  conduction  into  the  gas.  As  the  frequency  increases,  more  heat 
energy  is  conducted  from  the  particle  and  the  temperature  variation  at  the  surface  decreases  and 
changes  phase  with  respect  to  the  laser  excitation.  Chan12  has  discussed  this  effect  in  regard  to  thermal 
blooming  from  high  intensity  laser  beams  in  aerosols.  When  this  term  is  neglected,  the  ratio  of  the 
heat  and  mass  term  in  the  source  equation  is  frequency  independent.  The  inclusion  of  Eq.  (60c)  in  the 
formulation  provides  a  diagnostic  for  experimental  investigation  of  mass  transfer  on  photoacoustics. 
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Figure  4a.  Magnitude  of  the  ratio  f0{ji^rno)/  as 

a  function  of  neighborhood  radius  for  different  operating  fre¬ 
quencies  and  drop  sizes.  2000  Hz  and  10  micron  sized  drop 
(solid  line);  200  Hz  and  0.1  micron  sized  drop  (dashed  line); 
20  Hz  and  1  micron  sized  drop  (dotted  line).  The  ambient 
air  pressure  is  1  atm  and  the  ambient  temperature  is  290  K. 


Figure  4b.  Phase  of  the  ratio  /j(\/t2r?70)  as  a 

function  of  neighborhood  radius  for  different  operating 
frequencies  and  drop  sizes.  2000  Hz  and  10  micron 
sized  drop  (solid  line);  200  Hz  and  0.1  micron  sized  drop 
(dashed  line);  20  Hz  and  1  micron  sized  drop  (dotted 
line).  The  ambient  air  pressure  is  1  atm  and  the  ambient 
temperature  is  290  K. 


Figure  5  shows  the  magnitude  and  phase  of  the  acoustic  source  strength  for  1  pm  water  drops  in 
saturated  air  at  290  K.  In  this  figure,  we  set  Q0- 1 .  Also  shown  is  the  magnitude  and  phase  of  the 

photoacoustic  source  strength  per  unit  volume  for  an  identical  particle  with  no  latent  heat  (a  non¬ 
volatile  particle)  and  the  magnitude  of  the  individual  contribution  to  the  volatile  source. 


Figure  5a.  Magnitude  of  the  acoustic  source  strength  for 
a  single  drop  as  a  function  frequency  for  1  micron  sized 
drops.  Total  effect  (solid  line);  Heat  contribution  (dotted 
line);  Mass  contribution  (dash-dot  line);  Total  effect 
without  latent  heat  (dashed  line).  The  ambient  air 
pressure  is  1  atm  and  the  ambient  temperature  is  290  K. 


Figure  5b.  Phase  of  the  acoustic  source  strength  for  a 
single  drop  as  a  function  frequency  for  1  micron  sized 
drops.  Total  effect  (solid  line);  Heat  contribution  (dotted 
line);  Mass  contribution  (dash-dot  line);  Total  effect 
without  latent  heat  (dashed  line).  The  ambient  air 
pressure  is  1  atm  and  the  ambient  temperature  is  290  K. 
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At  lower  frequencies,  the  second  term  in  S  is  negligible  and  the  volatile  particle  source  strength 
is  smaller  than  the  non-volatile  source  strength.  The  cooling  of  the  particle  by  evaporation  reduces  the 
temperature  variation  at  the  drop  surface  and  therefore  reduces  both  the  amount  of  heat  conducted  into 
the  gas  and  the  amount  of  mass  injected.  At  high  frequencies,  the  cooling  of  the  drop  surface  due  to 
conduction  into  the  gas  becomes  dominant  over  cooling  due  to  evaporation.  The  magnitude  of  the 
thermal  source  alone  approaches  that  of  the  thermal  source  of  the  non-volatile  particle.  The  mass 
injected  into  the  gas  then  adds  into  the  heat  flux  term  to  increase  the  photoacoustic  signal. 

A  transition  from  the  source  strength  being  linearly  dependent  on  frequency  to  near  constant  as  o) 
is  increased  is  observed.  The  position  of  the  transition  depends  on  the  size  of  the  terms  in  Eq.  (55). 
The  second  term  is  the  only  term  that  is  sensitive  to  the  frequency  and  to  the  particle  size,  other  than 

the  leading  (r0/a)3  term  in  Eqs.  (57)  and  (58)  that  just  expresses  the  volumetric  fraction  occupied  by 

light  absorbing  material.  The  phase  of  the  photoacoustic  signal  varies  from  3tt/2  to  2n  as  the  second 
term  in  Eq.  (55)  increases  in  size. 


2.6  Applications 

These  equations  can  be  applied  to  estimate  the  effect  of  thin  layers  of  adsorbed  volatile  liquids  on 
the  photoacoustics  of  solid  particles.  It  is  expected,  for  example,  that  water  will  be  adsorbed  onto 
carbon  particles  at  vapor  pressures  below  saturation.  If  we  assume  fractional  coverage,  f,  of  the  solid 
by  the  liquid  and  the  heat  capacity  and  mass  of  the  liquid  is  small  with  respect  to  those  of  the  drop,  the 
photoacoustics  source  strength  can  be  determined  by  Eqs.  (57),  (58)  and  (59)  with  the  following 
substitution:  cpL  ->  cpS ,  the  specific  heat  of  the  solid,  and  pL->  ps,  the  density  of  the  solid.  In 

addition,  the  factor  f  is  introduced  into  Eq.  (49).  This  introduces  the  factor  f  as  reducing M  and  in 
reducing  the  third  factor  in  Eq.  (55).  This  also  assumes  that  the  solid  particles  are  spherical. 
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2.7  Conclusions 


A  theory  of  photoacoustics  for  volatile  aerosols  for  arbitrary  acoustic  frequencies  has  been 
developed.  Correction  terms  for  larger  particles  and  for  high  particle  densities  have  been  included  in 
the  analysis.  At  high  frequencies,  a  term  not  included  in  Baker’s  analysis  is  significant. 

The  contribution  of  evaporation  to  the  photoacoustic  measurements  of  light  adsorption  can  be 
determined  using  Eqs.  (55),  (57),  (58)  and  (59).  If  the  changes  necessary  to  describe  fractional  liquid 
coverage  on  a  solid  substrate  particle  are  incorporated,  the  source  strength  contains  four  unknowns, 
assuming  that  the  identity  and  properties  of  the  liquid-vapor  species  is  known.  Two  measurements  of 
magnitude  and  phase  of  the  source  strength  at  widely  spaced  frequencies  are  sufficient  to  determine 
whether  evaporation-condensation  is  occurring  on  the  absorbing  particles  and  to  evaluate  the  size  of 
that  contribution.  Multifrequency  measurements  would,  of  course,  provide  more  information. 
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